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Spinons are emergent fractional (spin-1/2) excitations of quantum magnets. While they are well 
established and understood in one dimension, characterizing them in two-dimensions, in theoretical 
models and real materials, has been challenging. Here we use Monte Carlo methods to study 
spinons in a resonating valence-bond (RVB) spin-liquid, as well as in a model which can be tuned 
from an antiferromagnet to a valence-bond solid (VBS); a crystal of singlets. We use a generalized 
valence-bond basis to compute the intrinsic size A of an individual spinon and the size A of a bound 
state. We confirm that spinons are deconfined, A — > oo and A finite, in the RVB. In the VBS 
A > A, i.e., the spinon is soft and shrinks as the bound state is formed. Both A and A diverge at 
the VBS-antiferromagnetic critical point. We conclude that the spinon deconfinement is marginal 
in the lowest-energy state in the spin-1 sector, due to weak attractive spinon-spinon interactions. 
Deconfinement should occur at higher energies. 

PACS numbers: 75.10.Kt, 75.10.Jm, 75.40.Mg, 75.10.Pq 



The concept of confinement originated in particle 
physics, where quarks are bound into hadrons by "gluon 
strings" , leading to a distance-independent force and the 
inability to isolate individual quarks 1 . This picture 
has been carried over to condensed matter physics as 
well. One example is the non-magnetic valence-bond- 
solid (VBS) state of a two-dimensional (2D) quantum 
spin system, whose ground state consists of a crystaline 
arrangement of two-spin singlets (valence bonds) . An el- 
ementary excitation corresponds to promoting one bond 
into a triplet. The members of this triplet (the spinons) 
stay bound to each other, due to a string of misaligned 
singlets forming when separating the spinons [2J. In 
a one-dimensional (ID) VBS the confining string does 
not exist, as there is no medium surrounding the phase- 
shifted dimers between the spinons. The spinons then de- 
confine and become independent elementary excitations. 
The spinon concept was originally introduced as a spin 
analog [3: of the fractional excitations of conjugated poly- 
mers [IJ, which are phase- twist solitons — domain walls 
between two out-of-phase dimerized chain segment. 

An ongoing quest in quantum magnetism is to identify 
2D systems in which spinons can deconfine. Deconfine- 
ment and weak confinement of spinons have been exper- 
imentally observed in one dimensional (ID) spin chains 
[5J and ladders [6], respectively, by neutron scattering. 
Deconfinement is characterized by a broad continuum of 
spin-1 excitations. These observations and conclusions 
are aided by the known triplet spectrum [7] of the Heisen- 
berg chain, where the existence of spinons is understood 
rigorously [8]. For 2D systems, this way of detection 
is difficult, due to a lack of established model spectral 
functions for spinons and the existence of significant con- 
tinua even in the magnon spectrum of antiferromagnets 
[5]. While alternative experimental methods have also 
been proposed [ID] , this remains a challenging issue. In 
numerical model studies, correlations of excess magneti- 
zation with impurities [111 112] and other inhomogeneities 



[13] have been used to study spinons. 

Here we use a quantum Monte Carlo (QMC) approach 
for homogeneous systems Q3] and address the mecha- 
nism of 2D spinon deconfinement promoted by Senthil et 
al. [15H17j . In their picture, a spinon in a square-lattice 
VBS can be viewed as vortex at the nexus of four VBS 
domains, with the core carrying spin-1/2. These vor- 
tices bind with anti- vortices to form the gapped spin-1 
triplons. The binding potential weakens as the fluctu- 
ations of the VBS increase and a deconfined quantum- 
critical (DQC) point is approached which separates the 
VBS from the standard antiferromagnetic (Neel) state. 

Deconfined spinons are expected in spin liquids, which 
are non- magnetic states with no broken symmetries |18j . 
The DQC point is an isolated, gapless (algebraic) spin 
liquid. In the near-DQC VBS, the triplon is predicted to 
become large, exceeding the correlation length [TCI 117] . 
and may be regarded as loosely bound spinons. Sev- 
eral QMC studies of J-Q models, where the Heisenberg 
exchange is supplemented with higher-order interactions 
|19j . support the existence of an unusual critical point 
|19H23j and the critical thermodynamic properties are 
also consistent with spinons 24 . 

Objectives and Findings — We study properties of the 
DQC spinons and their near-DQC bound states in the 2D 
J-Q model. Computing the spinon size and confinement 
length, we demonstrate that a naive picture of a large 
bound state of two small objects fails — the spinons them- 
selves are of size comparable to the bound state and are 
"soft" , shrinking when forming bound states. This dif- 
fers from the prototypical short-bond resonating valence- 
bond (RVB) spin liquid [551 i which we also study as a 
point of reference. Here the deconfined spinons are small, 
with radius of a few lattice spacings. 

Methods — In Ref. [O] we showed how spinons can be 
characterized for bipartite systems in the valence-bond 
(singlet) basis extended to S = 1/2 and S = 1 by in- 
troducing one (on a lattice with an odd number of sites) 
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FIG. 1: (Color online) Transition graphs for 5 = (a), 
5 = 1/2 (b), and 5 = 1 (c). The bonds and larger circles 
indicate valence bonds and unpaired (f) spins, respectively, 
in ket (black) and bra (red) states. The statistics of separation 
of unpaired spins in (b) and (c) can be used to characterize 
the spatial properties of spinons and their bound states. 

or two (for an even number of sites) unpaired spins. An 
5 = state is a superposition of states \Vn) in which 
all N spins are paired up into N/2 singlets between sub- 
lattice A and B sites, while for 5 = 1/2 and 5=1 
the basis states are, respectively, |Vjv-i) ® | to) and 
iViv-a) ® | to) ® I tb), where a and b are sites on differ- 
ent sublattices [37J HH] • Superpositions can be generated 
probabilistically either in variational states (amplitude- 
product states [29] and their generalizations [30 ) or in 
QMC simulations of Hamiltonians. Examples of config- 
urations are shown in Fig. [T] We use the formulation of 
the projector QMC method discussed in Ref. [51] . 

The valence-bond basis being overcomplete and non- 
orthogonal, one should not consider the unpaired spins 
directly as spinons. However, as shown in Ref. [14], to 
which we refer for details, the statistics of the separa- 
tion of unpaired spins in transition graphs (see Fig. [I]) 
are related to wave-function overlaps and can be used 
to extract the intrinsic size of a single spinon as well 
as the size of an 5 = 1 bound state. We denote the 
probability distribution between t spins in the bra and 
ket state of an 5 = 1/2 configuration -PaaM and be- 
tween those on different sublattices in an 5 = 1 state by 
Pab(t). The former contains information on the spinon 
size and the latter on the bound state (or lack thereof). 
The same information is contained also in the difference 
C(r) = Cs(r) — Co(r) between the z-componcnt spin-spin 
correlation function in the 5 = 1/2, 1 and 5 = states 
|14j . The method of statistics of separations is easier, 
since C(r) is statistically noisy for large systems. The 
latter method can, however, be extended to non-bipartite 
systems and other computational methods. 

Models — We use Monte Carlo sampling [29] to gener- 
ate valence-bond configurations for the short-bond RVB, 
as done in two recent studies to characterize its spin and 
dimer correlations [52] [55J- We use the RVB as a ref- 
erence point for an established U(l) spin liquid. It is 
expected to have deconfined spinons but their spatial 
characteristics have not been considered before, to our 
knowledge. Our main interest is to characterize spinons 
at the Neel-VBS transition of the J-Q model. Here we 
consider a Q-term consisting of three singlet projectors 
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FIG. 2: (Color online) Spinon distributions and correlations 
along the diagonal lattice direction of the RVB spin liquid. 
Black and red circles show the single-spinon overlap Paa{t) 
in an 5 = 1/2 state on a 65 x 65 lattice. The exponential 
decay (fitted line) gives the spinon size A = 2.96(1) lattice 
spacings. Red squares show the two-spin distribution Pab(t) 
in the 5 = 1 state on a 64 x 64 lattice, fitted to ~ l/r a with 
a « 0.6. The green triangles show the difference C(r) between 
the spin correlations in the 5=1 and 5 = systems. This 
quantity exhibits both a spinon-size effect (exponential short- 
distance decay) and deconfinement (weak power-law decay at 
long distances). There is a phase shift at r ~ 9y2. 

Cij = 1/4 - S; • Sj. The Hamiltonian is [H] 

H = — J Cij — Q3 CijCklCmn, (1) 

(ij) (ijklmn) 

where (ij) and (ijklmn) denote, respectively, nearest 
neighbors and three nearest-neighbor pairs in horizon- 
tal or vertical columns on the square lattice. The critical 
coupling ratio separating the Neel and VBS phases is 
9c = (Qi/J)c ~ 1-50 [21], and for large g the VBS order 
is very robust. 

We also study the ID J-Q3 model including static 
dimerization, with different couplings on even and odd 
bonds; 

H = — (J\C2i,2i+l + J2C2i+l,2i+2) 

i— even 

(2) 

i 

When Ji = J2 the quasi-ordered to VBS transition 
(which is of the same kind as in the frustrated Heisen- 
berg chain [34]) is at Q 3 /J = 0.1645 [14]. In this study 
we use this critical ratio and tune p = Ji j J\. The explicit 
dimerization introduced by p ^ 1 leads to confinement 
|12) . while the spinons are deconfined when p = 1. We 
will demonstrate similarities with the 2D J-Q3 model. 

RVB Spin Liquid — RVB spin liquids have been consid- 
ered as promising candidates for explaining high temper- 
ature superconductivity in cuprates when doped [23 ISS] • 
It is therefore also interesting to examine in detail proper- 
ties of the insulating host system. Recently, it was found 
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FIG. 3: (Color online) Spinon distributions for the 2D J- 
Q3 model on L — 128 and 129 lattices at g — 3 (VBS) and 
g c — 1.5 (critical). The lines are exponential fits for the g = 3 
VBS, giving the singie-spinon size A = 6.4 extracted from 
Paa(t) in the S =1/2 state and the confinement length A = 
3.1 extracted from P|j(r) in the 3 = 1 state. Both lengths 
diverge as g ~ > <? c , where the distributions decay algebraically. 

that the simplest equal-amplitude short-bond RVB is a 
quantum-critical VBS, with exponentially decaying spin 
correlations but power-law dimer correlations [321 I33| . It 
has been long expected that the RVB hosts deconfined 
spinon excitations [15] . 

The parent Hamiltonian of the short-bond RVB was 
found recently [36 . Although the S = 1/2 and S = 1 
states we study here may not be its exact lowest states 
in these sectors, one can still expect them to be good 
variational states — the actual excitations of the Cano- 
Fendley Hamiltonian should be very similar. 

We characterize the spinons in Fig. |J The S = 1/2 
distribution Paa{t) demonstrates a well-defined intrinsic 
spinon wave packet, decaying as e~ r '* with the spinon 
size A = 2.96(1). The S = 1 distribution Pas (7) is 
peaked at short distances and appears to decay as r~ Q 
with a ~ 0.6. This implies marginal deconfinement due 
to weak attractive spinon-spinon interactions. In Fig. [2] 
we also show that the length-scales observed in Paa(t) 
and Pab (r) are manifested in the correlation function 
C(r) as well. Hence, the distributions do capture actual 
physical, basis-independent length-scales [13] . 

J-Q model — The Neel-VBS transition has been de- 
bated for years [2J. In 2004, Senthil et al. presented a 
scenario encompassing several earlier works [37] [3S] 
and further proposed a mechanism leading to a generic 
continuous transition [T51TP7] . This scenario is at odds 
with the "Landau Rule" according to which transitions 
between the two ordered states should be first-order. A 
key aspect of the theory is that both order parameters 
arise out of spinons, which condense in the Neel state 
and confine in the VBS (where valence bonds can be re- 
garded as tightly bound spinon pairs). Exactly at the 
DQC point separating the ordered states the spinons 
should deconfine. Although opposing views have been 
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FIG. 4: (Color online) Divergent behavior of the correlation 
length £ and the confinement length A as the coupling ratio 
approaches the critical point (g c = 1.50) of the J-Q3 model. 
The curves are fits to the power-laws discussed in the text. 

put forward [3UM41] , the generic continuous nature of 
the Neel-VBS transition has support in QMC studies 
of J-Q [HI EH E21 El] and other J2U] models, including a 
predicted emergent U(l) symmetry reflecting the gauge 
structure of the theory, where spinons interact with an 
U(l) gauge field, as described by the non-compact CP 1 
field theory. 

The DQC scenario motivates us to investigate spinons 
directly, by computing the spinon size A and the confine- 
ment length A in VBS states and approaching the critical 
point. As shown in Fig. [3] both Paa(t") and Pab (?") are 
exponentially decaying in the VBS phase, with A = 6.4(1) 
and A = 3.1(1) at g = 3. Surprisingly, the intrinsic size 
of a single spinon is, thus, much larger that the bound 
state of two spinons. We interpret this as a softness of the 
spinon, which causes it to shrink when subject to attrac- 
tive interactions from an anti-spinon. This should be a 
signature of the vortex-nature of the spinon, as the oppo- 
site circulations of the members of the pair should lead 
cancellations away from the double-vortex core. Such 
shrinkage of vortices can also occur under certain condi- 
tions in superconductors [H] H3] . 

Approaching the critical point both A and A diverge, 
and at the critical point power laws Paa{t) ~ l/r a and 
Pab{t) ~ l/r 13 obtain, with a « /3 = 0.3 ± 0.1. In 
the DQC theory it is predicted that the divergence of A 
should be faster than the correlation length £; A cx £ 1+k , 
with k > and less than the exponent governing the 
rate of divergence of the cross-over length-scale of the 
emergent U(l) symmetry (which should be the longest 
length-scale in the system close to the DQC point) [16] . 
We extract £ from spin-spin correlations. The length- 
scales £ and A are graphed versus the coupling ratio in 
Fig. [4] Using lattices sizes L up to 128, we can reliably 
extract A and £ when they are roughly less than 10 — 
beyond which the size-dependence becomes significant 
and extrapolating to infinite size is difficult. Although 
we can therefore not reach far into the asymptotic scaling 
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regime, the dependence on g — g c is still consistent with 
the expected power-law divergence, if we allow a constant 
correction, i.e., fitting to the forms A = a + b(g — g c )~^ 
and £ = a + b(g - g c )~ v ■ We then find p = 0.7(1) and 
v = 0.8(1). The correlation-length exponent v — 0.60(2) 
was obtained in Ref. [25] based on finite-size scaling col- 
lapse for larger systems in the close neighborhood of the 
critical point. The results based on Fig.[4]may be affected 
by further non-asymptotic corrections. Regardless of the 
precise values of p and i^, it is clear that k = fi/v — 1 is 
very close to 0. This is consistent with the value 0.2(1) 
obtained in Ref. [U] for the exponent governing the U(l) 
to Z4 cross-over. 

ID deconfinement — In the ID VBS phase, without en- 
forced dimerization, J\ = J2 in Eq. ([2|, the spinons are 
small and deconfined, as shown in Ref. |14j . By turning 
on a symmetry-breaking dimerization, p = J2/J1 > 1, 
one can tune the confinement length from arbitrarily 
large to arbitrarily small (as will be discussed elsewhere 
[4"4"]). Here, to compare with the 2D model approaching 
its critical point, we instead show results for g = Q3/J 
fixed at the critical value g c — 0.1645 when J\ = J2. 
Keeping g = g c and turning on the static dimerization, 
p > 1, we observe in Fig. [5] that A sa A (A being slightly 
larger), with both lengths diverging as p — > 1. This is 
similar to the behavior observed in the 2D model (apart 
from the lack of spinon shrinkage). Thus, the fact that 
the spinon size and the confinement length are both di- 
vergent does not invalidate the deconfinement picture. At 
the critical point, the ID "spinon shape" distribution 
Paa(t) also decays as a power-law, while the pair distri- 
bution Pab(^) is peaked at the longest distance, reflect- 
ing marginal (critical) spinons subject to weak repulsive 
interactions [14] . The main difference in the 2D model is 
that the effective spinon-spinon interactions are attrac- 
tive, not only in the VBS phase but also at criticality. 

Conclusions — We showed that a spinon in the 2D J- 
Q3 model shrinks when a bound-state (triplon) is formed. 
This is most likely a generic aspect of VBS states; a con- 
sequence of the vortex nature [T7] of the spinons. We 
should stress here that the reason why the vortices do not 
annihilate is that we restrict the system to the 5=1 sec- 
tor, while spinon annihilation would bring it back to the 
5 = ground state. Both the spinon size and the con- 
finement length diverge as the critical VBS-Neel point 
is approached, and at the critical point the distribution 
functions decay as power laws. This is also necessary for 
continuity, because in the Neel state both distributions 
become flat (as we discussed for ID systems in Ref. [14] 
and have also verified in 2D), when the spinons com- 
pletely loose their identity as individual objects. 

Our scenario deviates from a simple picture of the near- 
critical triplon being a large object formed by two small 
particles. The question then is: Are the spinons nev- 
ertheless deconfined, independently propagating excita- 
tions? We showed with a known example that in ID that 
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FIG. 5: (Color online) Spinon statistics of the J1-J2-Q3 chain 
with Q3/J1 = 01645 (the critical point when Ji = J2) for 
three ratios p = JijJ\. The straight lines for p > 1 are 
exponential fits and the curves for p = 1 are power-law fits. 



is possible. The most plausible scenario in 2D is that the 
behaviors found here are due to weak attractive interac- 
tions between a spinon and an anti-spinon — as would in 
fact be expected based on the DQC theory, where inter- 
actions are mediated by a gauge field (but the sign of the 
effective spinon-spinon interactions is not immediately 
clear). Note that we have only studied the lowest 5 = 1 
state, and in higher states the kinetic energy should over- 
come the weak attraction, leading to essentially indepen- 
dent spinons, as found when comparing the critical J-Q 
model with a gas of free bosonic 5 = 1/2 excitations |24j . 

In the simple RVB spin liquid, we also found weak 
attractive interactions, but in this case the spinons are 
small. The differences between the RVB liquid and the 
critical J-Q3 model should be due to the fact that the 
spin correlations decay exponentially in the former (and 
there is a spin gap), while they have a power-law decay 
the latter. Both models exhibit algebraic dimer correla- 
tions, which should be responsible for the residual spinon 
interactions. In the ID case, full deconfinement can be 
seen thanks to repulsive interactions, although the criti- 
cal spinons themselves are not small particles (as in the 
2D J-Qz model), being instead marginally localizable ob- 
jects described by power-law overlaps. Still, it is rigor- 
ously known that these marginal particles do propagate 
as individual degrees of freedom [S] . 

Our results suggest that the critical J-Q3 model and 
other related systems (generic Neel- VBS critical points) 
are close 2D analogues to the ID critical spin chains, with 
the differences essentially due to the different topological 
aspects of the spinons; vortices versus kink solitons. 
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